Light Transport in Random Media with "PT-symmetry 
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The scattering properties of randomly layered optical media with PT-symmetric index of refrac- 
tion are studied using the transfer-matrix method. We find that the transmitance decays expo- 
nentially as a function of the system size, with an enhanced rate ^(W) -1 = £,o{W)~ + £,-y{0)~ , 
where £o(W) is the localization length of the equivalent passive random medium and £ 7 (0) is the 
attenuation/amplification length of the corresponding perfect system with a PT-symmetric refrac- 
tion index profile. While transmitance processes are reciprocal to left and right incident waves, the 
reflectance is enhanced from one side and is inversely suppressed from the other, thus allowing such 
PT-symmetric random media to act as unidirectional coherent absorbers. 

PACS numbers: 42.25Dd, 11.30Er, 03.65.Nk 



Introduction-Wave propagation in naturally occurring 
or engineered complex media, is an interdisciplinary field 
of research that addresses systems as diverse as classical, 
quantum and atomic-matter waves. Despite this diver- 
sity, the wave nature of these systems provides a common 
framework for understanding their transport properties. 
One such characteristic is wave interference phenom- 
ena. Their existence results in a complete halt of wave 
propagation in random media which can be achieved by 
increasing the randomness of the medium. This phe- 
nomenon was predicted fifty years ago in the framework 
of quantum (electronic) waves by Anderson PQ and its ex- 
istence has been confirmed in recent years in experiments 
with matter [2] and classical waves [3HZ]- 

While the localization of classical waves has been well 
understood by now [5], only in the last decade has light 
propagation in active random media been pursued in- 
tensively [9 16J. Due to the absence of a conservation 
law for photons, light may be absorbed or amplified in 
the medium while phase coherence is preserved. This 
interplay of absorption or amplification and localization 
has been studied by using the Helmholtz equation with 
an imaginary dielectric constant of an appropriate sign. 
Several interesting results have been found, such as the 
dual symmetry of absorption and amplification for the 
average transmitance and the localization length [14l [15] , 
the sharpness of back scattering coherent peak and the 
statistics of super-reflectance and transmitance [9"HT3] . 

Quite recently, the possibility of synthesizing a new 
family of artificial optical materials that instead rely on 
balanced gain and loss regions has been suggested [T71 - 
[24] . This class of optical structures deliberately exploits 
notions of parity (P)and time (T) symmetry [2"5T - f2"T] as a 
means to attain altogether new functionalities and opti- 
cal characteristics |17) . In optics, VT symmetry demands 
that the complex refractive index obeys the condition 
n(r) = n*(—r), in other words the real part of the refrac- 
tive index should be an even function of position, whereas 
the imaginary part must be odd. VT symmetries are 
not only novel mathematical curiosities. In a series of 
recent experimental papers VT dynamics have been in- 



vestigated and key predictions have been confirmed and 
demonstrated [19] [21] [28] [29] . These include among oth- 
ers, power oscillations [IT] [TH] I23J 5 absorption enhanced 
transmition [21) . double refraction and non-reciprocity 
of light propagation |17j . In the nonlinear domain, such 
pseudo-Hermitian non-reciprocal effects can be used to 
realize a new generation of on-chip isolators and circu- 
lators [22 . Other results within the framework of VT- 
optics include the study of Bloch oscillations [30], the 
realization of coherent perfect laser absorbers [31] and 
nonlinear switching structures |32j . 

Work has also been done on disordered VT systems, 
with main focus on the spectral properties of the corre- 
sponding PT-Hamiltonians [33J 131] ■ However, relatively 
little has been done concerning their transport proper- 
ties. In this paper, we will examine the transmitance and 
reflectance through one-dimensional (ID) 'PT-symmetric 
systems with random index of refraction (see Fig. [T]) . We 
show that the exponential decay rate of transmitance, 
which defines the inverse localization length is asso- 
ciated with the harmonic sum of the localization length 
£o(W) of a passive system with the same degree of ran- 
domness and the amplification length £ 7 (0) of a periodic 
PT-symmetric system with the same degree of gain/loss. 
Furthermore, we find that the asymptotic value of the 
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FIG. 1: (Color online) A one-dimensional PT-symmetric mul- 
tilayered random medium. The gain/loss refraction index pro- 
file is uniform (see Eq. |TJ) with the loss side on the left (light 
color) and the gain side (dark color) on the right of the struc- 
ture. The real part of the refraction index contrast n R is ran- 
dom, uniformly distributed around no, and n R (— z) = n R (z). 
For large enough system sizes (or strong disorder and/or large 
gain/loss) the system acts as a high performance absorber if 
the incident wave is entering from the lossy side of the struc- 
ture (light arrows), while it "super-reflects" if the incident 
wave enters the structure from the gain side (dark arrows). 
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reflectance follows a single parameter scaling law which 
is dictated by the ratio A = £a(W)/£y(0). Finally, we 
show that while the transmition processes are reciprocal 
to left and right incident waves, the reflection is enhanced 
from one side and is inversely suppressed from the other, 
thus allowing such "PT-symmetric random media to act 
as unidirectional coherent absorbers (see Fig. [IJ. 

Mathematical Model-We consider a one-dimensional 
(ID) active disorder sample having a random VT- 
symmetric refractive index distribution n(z) = no + 
n R (z) + in l (z) in the interval \z\ < L/2. The system 
is embedded in a homogeneous medium having a uni- 
form refractive index n for \z\ > L/2 (see Fig. [I]). 
Without loss of generality, below, we will assume that 
the refraction index n outside the disordered medium 
is no = 1. Here n R represents the real index contrast 
and n 1 the gain/loss spatial profile. In experimental re- 
alization in optics [171 H9j , these amplitudes are small, 
e.g. n R ,n l -C no- For simplicity, we will assume that 
the sample is composed of an even number, L, of lay- 
ers of uniform width d, each with a constant real refrac- 



tive index n K (zj — d/2 < z < Zj + d/2) 



given by 



a random variable with a uniform distribution between 
(— W, + W) which satisfies the 'PT-symmetric constraint 
n R (zj) = n R (—Zj). Specifically we will assume 

n(z) = n + n R (z) + for - L/2 < z < 

= n + n R (z) - i 7 for < z < L/2 (1) 

where 7 > is a fixed gain/loss parameter n l . Although 
the majority of our simulations below have been done 
for n l (±z) = =F7, we have also checked that our results 
apply for the case that n I (±z) = =F7 + 5n / where Sn 1 is a 
random variable given by a uniform distribution centered 
at zero. Since the qualitative features remain the same, 
we will not distinguish between these two cases. In this 
arrangement, a time- harmonic electric field of frequency 
u! obeys the Helmholtz equation: 



d 2 E(z) 
dz 2 



+ ^n 2 (z)E(z) = 



Eq. (j2| admits the solutions E (z) — E f e + E b e 
for z < -L/2 and E+(z) = E^ 



(2) 



—ikz 



+ „ikz 



+ g—ikz 



for z > 



L/2 where the wave- vector k — nouj/c. The amplitudes 
of forward and backward propagating waves outside the 
disorder domain are related via the transfer matrix M: 



E 

E/ 



f 



Mn Mi 2 
M 2 i M 22 



e: 



(3) 



The transmission and reflection amplitudes for left (L) 
and right (R) incidence waves, can be obtained from the 
boundary conditions E^ = (EJ = 0) respectively, and 



are defined as th 



E f 



Bi 



These can be expressed in terms of the transfer matrix 
elements as t L = t R = t = ^\r L = -jj^^r = ^ 
[351136]. While the transmitance T = \t\ 2 for left or right 
incidence is the same, this is not necessarily the case for 
the left and right reflectance Rl = |^l| 2 and Rr = \rn\ 2 
respectively. Furthermore, from the above relations one 
can deduce that r L r* R = (1 - |i| 2 ) [3H [37]. Thus one 
establishes the conservation relation [31] [37] : 



\/ RlRr 



11 -n 



(4) 



When considering Hermitian scattering systems, where 
Rl = Rr and T < 1, this relation recovers the familiar 
T + R = 1. The non-Hermitian systems studied in this 
paper, however, exhibit non-identical reflectances, and 
Eq. Q describes the connection between the two. 

Below we investigate the scaling properties of transmi- 
tance T and left and right reflectances Rl , Rr from such 
optical structures, with respect to the disorder strength 
W, the gain/loss parameter 7, and the wave- vector of the 
incoming wave. We have used various random refraction 
index contrasts W € (0,0.5), and gain/loss parameters 
7 G (0, 0.01) which is typical for optical media. We use 
slabs with L = 10 to L = 10 4 number of layers, each 
having width d = 1. The logarithmic averages (InT) and 
(lni?) are performed over 10 4 disorder realizations. 

Transmittance-The transport properties of passive (no 
gain or loss) disorder systems have been thoroughly stud- 
ied. At large length-scale, such systems exhibit an expo- 
nential decay in transmittance (see middle black line in 
the inset of Fig. [2|. The associated inverse decay rate 
£o(W), reflects the degree of randomness and it is defined 
as 



V£o(W0 



lim (lnT)/L 



(5) 



For ID random media we have that £o(W) ~ 1/W 2 0. 

The other limiting case of an ordered 'P7'-symmctric 
medium, can also be treated analytically. One can explic- 
itly solve for the electric field inside the perfect VT lay- 
ered structure subject to scattering boundary conditions. 
The resulting expression for the transmitance reads: 



T = 



8(l + 7 2 ) 2 

7 2 (4 + 7 2 ) (j 2 cosh(2fcL7) - cos(2fc£)) + 8 + 57 4 - 7 6 + 4 7 2 (l + cos(fcL) cosh(fcL7) - 7 (2 + 7 2 ) s\n(kL) sinl^fci^)) 



(6) 



For large i-values, the term involving cosh(2fci7) be- comes dominant and the transmitance decays exponen- 
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tially as shown in the inset of Fig. [2] (see upper red line) . 
For the experimentally relevant case 7 <C 1, the asymp- 
totic decay in transmitance can be found from Eq. ([6]) 
to be 



16e -2fc 7 L 

7 4 (4 + 7 2 ) 



The corresponding decay rate of transmitance is 



1 



£r(0) 



= — lim 

L— >oc 



1 



In T -> 2fc 7 



(7) 



(8) 



which can serve as an operative definition of the so-called 
attenuation/amplification length £ 7 (0). 

On the other hand, for system sizes L smaller than a 
critical length-scale L c the transmitance remains approx- 
imately constant T«l. Near the critical length L sw L c , 
large oscillations in the transmitance emerge (see upper 
red line in the inset of Fig(2]) after which the transmi- 
tance decays according to the expression given by Eq. 
0. The value of L c can be evaluated approximately by 
the condition Too(L = L c ) = 1 which leads to the follow- 
ing expression 



1 



In 



( 16 ) 



2fc 7 ~ V 7 4 (4 + 7 2 )- 



(9) 



The existence of a critical length-scale L c is character- 
istic of gain media and is associated with the lasing 
threshold for which T diverges. Below this length stim- 
ulated emission enhances transmitance through the gain 
medium. On larger length scales stimulated emission re- 
duces transmitance. The slopes of InT at both sides of 
the maximum are approximately symmetric. In contrast, 
in the case of a PP-symmetric refraction index the in- 
crease of the transmitance for L < L c which is due to the 
gain, is balanced by the equal amount of loss which is 
symmetrically arranged inside the medium. As a result 
T « 1 for L < L c . Nevertheless, this gain/loss balance, is 
not able to smooth out the diverging behavior of T near 
the lasing threshold (see inset of Fig. J2j. 

Let us finally consider the case of PT-symmetric dis- 
ordered slab geometry. A representative behavior of the 
transmitance T as a function of the system size L is 
shown in the inset of Fig. [2] (see lower green line). To 
understand the exponential decay of T, one needs to 
consider the simplified geometry with index of refraction 
given by Eq. (JlJ. For the case of Eq. (JlJ the transfer 
matrix of the total PT-symmetric system is the prod- 
uct of the transfer matrix Mi associated with the lossy 
sub-system and M g associated with the gain sub-system. 
The corresponding transmitance through the combined 
system is given by 



T 



\T,T, 



i 1 g\ 



1 - nr q 



(10) 



It is thus sufficient to know the scaling behavior of each 
of the terms on the rhs of Eq. ( 10 ) in order to pre- 
dict the scaling behavior of T. These terms have been 



studied in Refs. [T5], where it was found that both ab- 
sorption and amplification lead to the same exponential 
decay of the transmitance which in both cases is en- 
hanced with respect to a passive disordered medium by 
the strength of the gain (or loss) rate i.e. (lnT;. g ) = 
— (2&7 + £o(W)~ 1 )L/2. Somewhat counter-intuitively, 
the sample with amplification also exhibits exponentially 
decaying transmitance due to the enhanced internal re- 
flections from the boundaries. Using the duality relation 
|15) r/r* = 1 for the reflection of an amplifying or at- 
tenuating medium (with the same rate of gain or loss 
respectively) applied for L/2 > £o{W) we get 

(In T) = (lnT l ) + (lnT g )-2(ln\l-r l r g \) ( n ) 
= (lnT,) + (lnT fl )-2<ln|2(l-cos(20))|) 

where 9 is the phase of the reflection amplitude r;. As- 
suming that 9 is a random variable uniformly distributed 
on the interval [0, 2ir] [IB], we get that the last term af- 
ter performing the average over the random variable 9 is 
finite. Therefore we get: 



lim 



(InT) 



==-(2k 1 + t (W)- 1 ) 



(12) 



From the above argument, we conclude that the localiza- 
tion length for a 'PT-symmetric disorder medium is: 



(13) 



In Fig. [2j we show the results of our numerical simula- 
tions for a disordered PT-symmetric sample. The ex- 
tracted localization length nicely follows the scaling be- 
havior indicated by our theoretical arguments. 

Reflectance- We proceed with the analysis of the re- 
flectances. In the case of random media with only gain 
or loss it was found in Ref. jT5] that the reflectances of 
two distinct disordered systems, one with gain strength 
—7 and the other with loss strength 7 satisfy the follow- 
ing reciprocal relation between them: 



RqainRloss 1 



(14) 



Specifically it was found that for gain media i? ga in > 1, 
while for lossy media we have via Eq. ( 14 ) the reciprocal 
behavior i?i oss = ^g a i n < ^ is important to stress here 
that such systems do not distinguish between left and 
right incidence, that is Rl = Rr for each of the cases. 

On the other hand, PT-symmetric systems distinguish 
the reflection between left or right incident wave, that 
is, Rl ^ Rr in general (see previous discussion). This 
phenomenon has already been observed for periodic VT- 
symmctric structures in Ref. |24l . Moreover, in the pres- 
ence of random index of refraction, we have previously 
concluded that the transmitance is effectively diminished 
exponentially with a rate 1/£ 7 (W) given by Eq. (13 1. 
Using the conservation relation Eq. Q we get 



RrRl —> 1. 



(15) 
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FIG. 2: (Color online) The numerically extracted localiza- 
tion length £ 7 (W) for various gain/loss parameter 7 (not in- 
dicated in the figure) is plotted, rescaled with £o(WO versus 
the scaling parameter ^o(W)/^(0). The symbols and colors 
indicate different wavelengths k of the incident wave, and dis- 
order strengths W. The black line indicates the theoretical 
prediction of Eq. ( 13 1. The meshed symbols correspond to 



some typical £,~/(W), for the scenario where the imaginary part 
of the refraction index is n 1 = 7 + Sn 1 where Sn 1 is a random 
variable uniformly distributed around zero. In the inset we 
report (InT) against the system size for: a system with only 
gain/loss and n R = no constant (upper red); a random layer 
medium with n 1 = (middle black); and a PT-symmetric 
disorder medium (lower green). 



Although this relation is similar to Eq. (14), it should be 



emphasized once more that in the case of PT-symmetric 
disorder media the medium behaves simultaneously as 
a gain medium (i.e. having > 1) and as a lossy 
medium (i.e. it can enhance absorption of incoming co- 
herent waves R R < 1). The reciprocity of the left and 
right reflectances is clearly demonstrated for some repre- 
sentative cases in the inset of Fig. [3j 

A natural question is associated with the scaling be- 
havior of the asymptotic value of the reflectances Rf? R as 
a function of the disorder strength W and the gain/loss 
parameter 7. We speculate that a one-parameter scaling 
law describes the asymptotic reflectance i.e. 

RT R {l,W) = f{k), where A = £ (W)/£ 7 (0) (16) 



L R is the geometric mean of asymptotic re- 



Here, R^ j: 

flectance, that is, R^ R = exp((\nRf R )). We have tested 
our hypothesis numerically. To this end, we have ex- 
tracted R^ R from our data for various values of 7 and 
W and plot them against the scaling variable A. The 
results are presented in the main part of Fig. [3] We 
find that for realistic values of the gain/loss parameter 
7 < 10~ 2 — 10~ 3 the data nicely follow the one-parameter 



scaling hypothesis (16). As the scaling parameter A in- 



gain/loss parameter 7), the asymptotic value R^ R de- 
creases/increases. Such a behavior allows us to use the 
proposed structure as a unidirectional coherent absorber 
that can increase absorption by tuning up the scaling 
parameter A. We want to mark that our structure, is 
different from the one suggested in Ref. [35], where it 
is shown that a disordered system with a single absorb- 
ing element causes coherent enhanced absorption if the 
phases of the input field are appropriately manipulated. 
Instead, we are addressing a different problem, where we 
have broadband absorption in one direction without a 
need for phase manipulation of the incoming wave. 

Conclusions- We have investigated the transport prop- 
erties of one-dimensional 'PT-symmetric disordered lay- 
ers. We have found that the localization length 
£ 7 (W), defined as the inverse decay rate of the trans- 
mitance, is smaller than the localization length of the 
passive disordered system fo(W) and from the ab- 
sorption/amplification length £ 7 (0) of a periodic PT- 
symmetric medium. At the same time the reflectance 
depends on the direction of the incident wave: while for 
incident waves entering the medium from the gain side 
it is enhanced, it is suppressed if the wave enters the 
medium from the lossy side. The reduction/enhancement 
of the reflectance is dictated by a one parameter scaling 
A = Co(W / )/^ 7 (0) and allows us to use such structures as 
unidirectional quasi-perfect coherent absorbers. 
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FIG. 3: (Color online) Insert: Typical reflectances .Rl,r ver- 
sus the system size (number of layers L). The main figure 
displays the asymptotic value of | log R°° j versus the scaled 
parameter A. For large values of A the |logJ?°°| increases 
indicating that (i.e. reflectance for an incident wave en- 
tering the sample from the lossy side) diminishes. In this 
domain, the sample acts as a unidirectional absorber. We use 
the same symbol and color coding for our data as in Fig. [2] 
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